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ABSTRACT 


Evaporation from Lake Wabamun, Alberta is estimated over a 
Ponvods Of 5e/S Sdays during the summer vote 2/4eusing whourly 


data in a mass-transfer equation of the form: 


Fils p& Cu-u, MSY SD 


It is shown that a wind-and-stability-dependent formulation 
EOre sce provides an estimate for evaporation which does not 
differ significantly from evaporation calculated using a 
constant value for C, = 0.0016. The evaporation calculated 
using the mass-transfer method is found to be consistent 


With lake evaporation estimated by applying a Simplified 


water budget to the lake. 


The use of data in the form of daily averages of wind, 
temperatures, and specific humidities in the mass-transfer 
equation yielded a 10% underestimate of evaporation over the 


period of study. 


Upper and lower bounds to enhanced evaporation from 
Lake Wabamun due to the presence of heated plumes from power 
plants are formulated, and it is found that slightly less 


than one half of the excess thermal pollutant goes into 


TOAST RAE 


A) teva tb tegsshorE- FSISUEA . UROL ened ORs no. sezoqgeve 

7 

Fusimih oes py: 3¢ hotaeg) « 
7 

senas “st 30. Fess *ciedii?-sepr 6 Ge e7s6 ‘ 


7 


: is 
pol oe ha Le i ’o-geecoaiw 4-354) deude Bt SEY 
ae ni del ir ‘ ay » i & i ‘es Te eed a 7S *~ s 162. 

rebne fotttestiy ivkregaeh> eer SiiedateReete astaee 
: IMAI 4 ee Sd = Le ay -aeta i | 7 Pe ae Pp ~wiiay ea8 HOR -_ 


See Amy - vel m9 rane at AIS Yodehaly -Beau ai? isga 
4 
bin y>4 fen piv }? wil & Sx PLIES GO Pay “AS 4 ve wsad 


—rtel wie vt Spee 710d" 


(Bale 76 wepézeva G1 P64 3o's10' aid pl pied 1° Oey eRe, 
Wehiieti-elbw afk Ad EWDPeNWed vi Blow Ane  eetwese 
ad Teva ADS TEAUTS Bi ty WIREES SHI Gi cor 4 cape 8 ; 
| oo. gute me he 


vy 


Daas: > “i iow 
' _ ph! sg ii zy 7 fara 


a saps = rots 


z 
~~. 


- 2m | 7 He mn 7 * A incoent * 
| hel) Pa ed ate 
‘ Pp 


aa} 


evaporating water, while the remaining energy is er ered 
between heat loss due to radiation and sensible heat 
transfer. It is estimated that the power plants operating 
at a capacity of 900 MW would increase annual evaporation 
over the lake by approximately 2 cm. The enhanced 
evaporation is found to be less than the computed mean 


annual discharge (6 cm) from Wabamun Creek. 
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CHAPTER TE 
INTRODUCTION 
1.1 Background for the Study 
Lake ¥abamun, Alberta has been the site # cof 


environmental impact studies in recent years. The lake is 
Within an hour's drive of Edmonton,-and serves that centre 
and many towns and districts as a recreation area {Fig 1). 
In addition to summer cottages, beaches, the Paul Band 
reserve, and a provincial park bordering its shores, there 
are two power plants located at the eastern end of the lake. 
The power stations use the lake water for coolant, and 
inject heated plumes of water onto the lake surface. These 
plumes have been a concern of the environmental sciences for 
a number of years. Along with claims of accelerated weed 
growth and fish kills, voiced by local residents, there has 
also been some Suspicion that, in time, the increased 
evaporation resulting from the heated water would lower 


water levels in the lake. 


Hogg Cho) using fast-response instrumentation 
compared evaporation over the hot plume, and over the 
unaffected lake surface, found greatly enhanced evaporation 


over the heated water. Hogg could not, however, make firm 
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conclusions about the long-term effects, Since his findings 
were based on two warm days in August and a cold day in 


Ockober. 


in response to public concern, Calgary Power Limited 
has constructed a large cooling pond and has diverted North 
Saskatchewan River water to service the Sundance plant on 
the south shore of the lake. In addition, water levels at 
Lake Wabamun in 1974, as in many Alberta lakes that year, 
were extraordinarily high and public concern for enhanced 


evaporation has waned somewhat. 


Thevstudies of Fritz and Krouse — (19/3), and Nursall, 
Hiwatt Zee OO Lec hee usOLOpICceCONStLLUeNtS and 
salt content of Lake Wabamun indicate nha the lake may have 
a Significant throughflow and inflow of groundwater. Their 
hypothesis is supported by Carlson (1970) who suggests that 
a buried Pleistocene riverbed connects Lake Wabamun with 
Lake Isle, and drains some of this water to the south east. 
TheLte thesis SisSee that. thus aquitern both charges, and 
dischatges, Water trom the lake. They Trutther estimate that 
thebe 2S a netwinflow of groundwater into the lake which 


would be sufficient to raise the water level of the lake by 


3-6 cm monthly in the absence of evaporation and rainfall 


It was the hope of this study that, with data collected 
during 1973 and 1974 from a meteorological station network 
on the shores of Lake Wabamun, the effects of the heated 


plumes on evaporation, and the net flow of ground water into 
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the lake could be calculated. No comparable studies 
involving meteorological JakeShore measurements have been 
made for any Alberta lake. All prior estimates for lake 
evaporation in Alberta were based on data from existing 


climatological and synoptic observing networks. 


eS —— 


Many methods are available for the study of mass and 
heat transfers . A battery of fast-response instrumentation 
could provide a good and reasonably representative estimate 
of lake evaporation by eddy-correlation methods, but the 
resources of equipment, and the labour required for such a 
study are prohibitive. Atmometers such as evaporation pans 
serve as analogues to lake conditions, and through the vse 
OfrmLOLriwlas SUCH as those Ot.) KOhISr et val. (1955) and 
others produce reasonable estimates of lake evaporation from 
Shallow lakes and reservoirs. The mass-transfer approach 
for estimation of evaporation (Jacobs 1942 and others) was 
used in this study, because it provides a direct check of 
Hage's (1974) calculation of Lake Wabamun evaporation based 
on monthly averages of meteorological variables. The mass- 
transfer approach has an advantageously low cost in terms of 


labour and instrumentation compared to other methods. 


Intercomparisons of lake evaporation calculated by 


differing methods were made by McKay and Stichling (1961). 
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Comparisons were made of evaporation estimated by the 
formulas of Meyer (1942), Penman (1948), Marciano and 
Harbeck(1954), and Kuzmin(1951) with evaporation estimates 
based on evaporation-pan, floating-pan evaporation, and 
water-budget measurements £OG YVeyburn Reservoir, 
Saskatchewan. Large discrepancies were observed between the 
various estimates of evaporation. Lteivass | pointed.out, 
however, by Harbeck(1961) that much of the observed 
discrepancy could have arisen from the manner in which the 
parameters of wind and water temperature were obtained for 
ERTSee Study « Buckler (1973) estimated evaporation rates for 
Lake Diefenbaker, Saskatchewan, a large deep lake, and 
concluded that stored heat played an important role in the 
evaporation process in the late fall. Spring and Schaefer 
(1974) used a modified form of the Lake Hefner mass-transfer 
eguation to obtain daily estimates of evaporation for Pabine 
Lake in Northern, BritisheColumbiate) Recently,., McBeans (19/5) 
computed fluxes of heat and moisture from Lake Ontario 
during a cold frontal passage using a mass-transfer equation 
to compute evaporation. A significant feature cf McBean's 
calculations was the use of a builk-transfer coefficient 


which varied with the wind speed. 


The Monin-Obukhov Similarity Theory (Lumley and 
Panotsky 1964) implies that the non-dimensional bulk 
transfer coefficients for momentum,heat and water vapour are 
Squaiweeco scach other) atthe limi or neutralestabi lacy. 


McBean (1975) used the results of Smith and Banke's (1975) 
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estimates of (C,) as estimates for C, to compute evaporation 
from Lake Ontario. Hicks (1972) evaluated the bulk-transfer 
coefficients for momentum, water vapour and sensible heat 
over water bodies of different sizes and found that, while 
drag coefficients increase slightly with wind speed at 
moderate to high winds, at low wind speeds they do not 
differ significantly from those over aerodynamically smooth 
surfaces. HiCks, salso found “ane increase Pots) the drag 


coefficient with increasing instability. 


After considering the constraints of available data, it 
was decided that a bulk-transfer approach would be used 
along with water-budget estimates of evaporation in an 
attempt to determine whether or not there exists a net 
groundwater inflow into the lake, and to estimate the 


effects of the thermal effluent from the two power stations. 


The formulation of the bulk-transfer coefficient 
(Jacobs 1942) used by Hage (1974) to estimate evaporation 
from Lake Wabamun made no ailowances for variations with 
Windspeed or stability of this coefficient. ain this thesis 
the theoretical results of Deardorff's (1968) calculation of 
grag coctricient variation @with Stabilitysiseused to ~adjust 
the transfer coefficient for varying stabilities over the 
lake. The wind-speed-dependent formulation for the neutral 
stability drag coefficient used in this thesis was based on 
Charnock's (1954) dimensional analysis argument, and the 


ebservations of Smith and Banke (1975), and Wu (1969) for 
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high wind speeds. Nikuradse's (1933) profile formulation 
specifies a bulk-transfer coefficient corresponding to 
aerodynamic smoothness for low wind Speeds. This 
formulation for the drag coefficient is in agreement with 
some of the results of Hicks (1972) . Chapter! is (oD) this 


thesis expands on the reasons for these choices. 


1.3 The Experiment 


— =e oe ae ee See 


There are few experiments which are carried out under 
idealiconditions, and this is not an exception. At the 


outset it was hoped to include data collected in the summer 


of 1973, as well as 1974 in the obServations. However, 
reduction Of “the. ‘wind “obServatirons for 1973 was mot 
COMDLEteCd An tine foe senciusion san thas “analysis. The 


conclusions found here were hased on 78 days. of data 
collected during the period 11 June 1974 to 28 August 1974. 
While the data were restricted to the summer months, some 
eStimates were made of the effects of the thermal plumes 
during the winter months on lake evaporation, and an 
estimate was made of the magnitude of the yearly evaporation 


from the plume using the bulk-transfer method. 
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DATA COLLECTION 


Lake Wabamun (Fig. 1)** Located? 767) “km west of 
Edmonton,is a large eutrophic lake surrounded by gently 
rolling parkland. The longest axis of this lake extends 
from WNW to ESE ‘more or less along the directions of the 
prevailing winds. The morphometry of Lake Wabamun is 


outlined in Table 1. 


The power stations and the observing sites are located 
at the eastern end of the lake (Fig 2). A thermohygrograph 
in a Stevenson screen, a Ryan water-temperature thermograph, 
a cup anemometer, a Stevens water-level recorder, and a 
staff gauge were installed on the north side of the lake, at 
the end of the Alberta Concrete Products Limited berm (Table 
2 and Fig 3). A thermohygrograph in a Stevenson screen, two 
cup anemometers, a Stevens water-level recorder, and a staff 
gauge were installed at the end of the intake canal for the 
Sundance Power Plant, on the south Side of the lake (Table 2 


and Fig 4). 


The topography of the shores adjacent to the instrument 
sites is shown in Fig 2. The steep slopes to the east of 
the Alberta Concrete Products Ltd. berm dictate that only 


the westerly to southeasterly wind directions off the lake 
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TABLE 1 


Morphometry of Lake Wabamun 
(aABLer ogg 1973) 
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Elevation APTA Gs) 

Area 82.5 kme 
Yolume ; | ca. 0.455 km? 
Length SIS PR. a 
Maximum breadth ; 626° Km 

Mean breadth 4.3 km 
Maximum depth as Wy 
Shoreline length Sie ek 


Area of surface drainage 372.4) kn < 
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Sundance Power Plant: 


FIGURE 2. 


Lake Wabamun ObServing 


Sites and Topography. 
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TABION 2 
Instrumentation 
Location Element Height Sensor 
(28/8/74) 
Sundance Wind speed 3.67 m Cup anemometer 
Wind direction Sis (on! Ga Vane 
Air temperature ero a Thermohygrograph 
Relative humidity 1.54 nm Thermohygrograph 
waterlevel Stevens recorder 
Wabamun Wind speed Sips We? air Cup anemometer 
Wind direction Sie je Vane 
Air temperature We. 0) 1 Thermohygrograph 
Relative humidity 1.05 mn Thermohygrograph 
waterlevel Stevens recorder 


Water temperature -6 cm Ryan Thermograph 
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FLGURE 3s 


Pictures of the Wabamun Observing Site. 
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FIGURE 4. 


Pictures of the Sundance Observing Site. 
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could be used for mass-transfer calculations. The 
topography immediately around the Sundance site was fairly 


flat, however. 


a i a ee ce ee ee ee ee eee 


The instruments (Table 2, Figs 3 and 4) were inspected, 
and calibration measurements were taken at approximately 
weekly intervals during the summer of 1974 by Mr.- Bob 
Weatherburn and Dr. John Honsaker of the Meteorology 
Divison at the University of Alberta. The cup anemometers 
were calibrated in the Department of Civil Engineering wind 
tunnel prior to, and after the experiments. The procedures 
followed for each instrument are outlined below. Because of 
the magnitude of error involved, the Assman psychrometer 


used for field calibrations also required recalibration. 


— ae eS 


AS a calibrating instuument, the Weather Measure Corp. 
Assman psychrometer model H 331 was a disappointment. When 
a@eecheck, twas ~madetin thes Laboratory of»the,yreadings of the 
wet and dry bulb thermometers used for calibrating the 
thermohygrographs, it was found that a discrepancy of 0.8 C 
existed between the wet and dry-bulb thermometers. The 
calibration of the Assman psychrometer was perrormed by 
placing both the wet-and dry-bulb thermometers initially in 


an ice bath and then, after the ice melted, observing the 
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temperatures as the water temperature rose to room 
temperature. BYSeSuUNtract Ing e906 CoPrLrom ac) ieeenetany— ou LD 
thermometer readings used in calibration of the 
thermohygrographs, the standard deviation between the wet 


ana dry—bul be thermometers was Leduced to 0.2) Ce 


2.2.2 Thermohygrographs 


The thermohygrographs were initially adjusted, and then 
cycled through changes in temperature and relative humidity 
within the environmental chambers of the Department of 
Botany of the University of Alberta. The calibrations 
applied to these instruments were obtained from weekly 


comparisons with the Assman psychrometer. 


Least-squares regressions of corrections of the forn: 


iF 2 Cet: + £60 
trae measured 17 (2. 1) 
R i yp laceek® as 
Micke, maszured 9 > CsA 
were made to the observations recorded by the 
thermohygrographs. The corrected relative humidity and 


temperature were plotted against time, relative humidity, 
and Pare ere to check for any remaining trends. No 
obvious trends existed, and there appeared to be no 
deterioration or aging of the hygrograph elements during the 


period in which the experiment was conducted. The 
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thermohygrographs were mounted Within an unaspirated 
stevenson screen, and thus there is the possibility of a 
considerable wind-speed-dependent time lag between the 
conditions in the environment, and screen measurements. 
Bryant (1968) examined the time-lag for the Stevenson screen 
and found lag times to vary from 30 min in calm conditions 
to 6.5 min with 7 m/s winds. Calibration measurements taken 
With the Assman psychrometer would be the instantaneous 
conditions of the environment. No compensation was ave for 


ee Os 


2.2.3 Ryan Water Temperature Thermograph 


The Ryan water-temperature thermograph located at the 
Wabamun Site was calibrated against periodic checks of water 
temperature using a standard thermometer. After addition of 
1-7 C, the Ryan temperatures agreed to within +0.9 C with 


the calibrating thermometers. 


TheLceawwelLeorNOAMOUDMZOUST I INStancesrs Whole ates coulda. be 
definitely concluded that the hot-water plume from the 


Wabamun power plant was influencing the water temperature 


readings. 


During the 78 days of the experiment water-temperature 
data were unavailable during the periods 0000 MDT, 29 June 


to 1700. MDT. 2 July. and 0600 MDT,» 5,August to, 1800 MDT.) 6 


Aug. 
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2.2.4 Stevens Water-level Recorder and Staff Gauges 


—we —. ee ee 


EXCeptsLOns Short gaps caused by jammed" charts etc. At 
the Sundance site, a nearly continuous record of lake levels 
waS available both at Wabamun and Sundance instrument Sites. 
The charts in the Stevens recorders were changed, and a 
reading was taken of the water level at each site from the 
staff-gauge at weekly intervals. Calibrations factors based 
on the site staff-gauge measurements were then applied to 
data extracted from each site. The calibration calculations 
resulted ina record of the height above water of the 
thermohygrograph units, used later in this work for mass- 


transfer calculations. 


The water-level measurements of the two Stevens 
recorders provide a “con¥annouss*record®%ef water-level 
observations. These records were combined to produce lake 
levels above mean sea level uSing water-level measurements 
Magssby-eNuttal (1974)-. Lt was also necessary toy ‘compensate 


for. wind setup in the water-level readings. 


A continuous record of wind speed and direction was 
available from the strip-chart record of the Rustrak event 
recorders. Readings from the strip-chart record were 
tabulated on magnetic tape, corrected for time and reduced 
to hourly means of calibrated” wind speeds using’ a 


minicomputer. The anemometers and their recorders were 
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calibrated in the Department of Mechanical Engineering wind 


tunnel using the procedures described by Honsaker (1975). 


2.3 The Local Weather Conditions During the Experiment 


—— eee es ES eee 


The unusually heavy snows of the winter of 1973-19741 
helped bo (produces «extraordinarily. «high! watepyelevels in 
nearby Lakes and sloughs during the spring and summer of 
1974. Water levels were high enough to produce a =Ilarge 


surface discharge from Wabamun Creek, at the east end of 


Lake Wabamun. 


Bhe summeurotmiosit, dopiae least thatapomerongupthin the 
18 \daysvot themstudy,ecouldgbepdescurbed £obethey masts part 
as ~Showery, and “cool. Guy 325 BOG tie 78 days measurable 
precipitation was recorded by the rain gauges. Except for 
localized or short-duration winds during storms, the monthly 
on winds were lower than normal. Relevant climatoiogical 
statistics2 for June, July, and August 1974 are given in the 
monthly summaries from the Atmospheric Environment Service 


below: 


June 


The month of June in Edmonton was mainly 
sunny, warm and contrary to month end Statistics, 
was generally dry. 


ae ee ee eee 


t Snowfall recorded at the Edmonton International Airport 


Was) 219. S.etuvsemia2 Gi nobmad Ly ; 
2 Since the summaries were written before metrication the 


units are non-metric. 
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Sunshine, wats 30 5eecCecorded) shouus, sawas, 040 
hours above normal and the sunniest since 1969. 


The mean daily temperature of 61.7 degrees 
was 3 degrees above normal. Temperatures tcranged 
from a low of 44 recorded on several days toa 
high of 88 degrees on the 23rd. 


LOtal proecipitataon of,458/ <inchesgeton ‘the 
month was much above normal. Of this amount, 3.15 
anches, or) 80.245 einches above the normal total of 
Z2uSteinches tora the entire month fell oneJdune,. 26. 
NEVO GOV en Ieee ndavdenour Traintalleancuntcrtor vuune 
were established on this date. 


Thunder recorded on 10 days, was above the 
normal of 4.4 thunderstorm days for the month. 


July 


Temperatures during the month of July in 
Edmonton were Slightly below normal during the 
first half but recovered to near or slightly above 
normal during the later half. 


Precipitation, for the gmonth totaled 4.97 
inches, 1.69 17ches above normal... Of thase “total, 
e226 inches febly during 10 of the first 12) days of 
the month. 


Two Significant storms passed through the 


area. A rain storm during the (eon 2 
accompanied by winds gusting to 51 mph dumped 2.24 
inches of Caine OveD sthew city scesuitange in 


Localizgedstloodings Locethe 2nd 2 memi ne lGrdays on A 
vicious thunderstorm crossed the northern portion 
of the city during the late afternoon of the 30th. 
Vindserat. “Chemindtstiia lead Doone DeaChedmgquste.or 
“7 mph however, near the storm's center, Namao 
reported gusts to 87 mph. Apart from these two 
windstorms average wind speed for the month was 
only 7.8 ° mph, the previous low average Speed Set 
abigh Pc Borers 


August 


The month of August in Edmonton was cool, dry 
and cloudier than usual. With a mean daily 
temperature of 58.5 degrees, = degrees below 
normal, the month was the coolest August since 
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1S OG. Temperatures reached into the 80's on only 
2 days. Those were on August 3 and 4 when the 
daily maximum temperatures were recorded as 85 and 
86 degrees respectively. The lowest temperature 
for the month was 42 degrees on the 13th. 


Precipitation for the month was 1.675 inches 
below normal and the first month Since July 1973 
that total precipitation was not above normal. 


Sunshine with 252.7 hours recorded was 16 
hours below the long term normal. Winds during 
the month were lighter than average. Highest gust 
was from the northwest at 45 mph on the 24th. 
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COAP TE Re Laer 


THEORY 


Most profile and mass-transfer theories used” in 
micrometeorology are derived by considering an infinite 
homogeneous level surface. For a lake the size of Wabamun, 
sucht liang appuoach jas justi fiedce Becauselthe Wakevas snail 
compared to synoptic-scale systems and because the layer of 
interest is Shallow, pressure-gradient and Coriolis forces 
may be neglected in the equations of motion. Since 
turbulent diffusion is much greater than molecular diffusion 
in the free atmosphere, molecular diffusion is neglected 
also. Furthermore, radiative flux divergence will be small 
except in cases of extremely stable stratification (Lumley 
and Panofsky, 1964). What remains is a constant-flux layer 
near the ground. With few exceptions, the stability regime 
over the lake favours lapse conditions The assumption of a 
constant-flux layer may not be fully justified over the hot- 
water plume, where horizontal gradients of water~-surface 


temperature exist. 


in summary, it is assumed that the earth"sS rotation, 


radiative flux divergence, and molecular diffusion, are 
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negligible. With the futher assumptions of steady state, 
horizontal homogeneity, and the Boussinesq approximation 
(density Saeecunct1on) sof | geopotential). onewsobteains the 


following equation of motion for the mean flow: 


> —p 
Ca naan cl ings BPS TEA (aire) 
thas = 
al ae e d2 
3.1.2 The Monin-Obukhov Similarity Theory and More 
Assumptions 
Liem Kernel yw Ole Starcing) pornte. for the transfer 
equations used in this thesis is the Monin-Obukhov 


Similarity theory. By defining the Monin-Obukhov length: 


3 
L = ~ Uy, 


eC) Ka 7) 


(3-2) 
where: 
LL = tel +> fF fey AYE Coa 


One may then represent the profiles of wind, temperature and 


humidity as functions of Z/L; 
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Garon) ae Ci) 


The bulk-transfer coefficients for momentum, heat and 


water vapour are defined as: 


z 

Hy Bib mener (nets (3.6) 
H 

Cece = ane ea Gees 
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é 
=i 
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= ed iN LG AGEN) (4-4, 

P 30) 
In’ order to remain consistent with the equations derived by 
Deardorff (1968) for the buik-transfer coefficients, a 
further assumption needs to be made. Deardorff has assumed 
that ces SANG de Keak - Thais 1S equivabens | to the 
assumption that ? = ¢, > ¢ lm = 9S, and that the 
roughness lengths for momentum, moisture and heat are equal. 
These assumptions are equivalent to considering the 
transport of momentum by pressure fluctuations to be 
negligible. Calder (1966) has reviewed the Monin-Obukhov 
Similarity theory, and che role of the pressure 
£luctuations. Businger et al. Gil) have obtained 
So =i eid. mk =0's3 0s These results would seem to 
indicate that this transport term should not be neglected. 


However, Dyer (1974) argues that the 10% reduction in wind 
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speed and 33% reduction in drag-plate shear-stress which 
Businger et al. Gigjj) used. torcalabrate theinmemstronents 
WasetoOO) aatge,se and that lower “values would result “in 
ae AAV aes : 

K,-1-9 with k=0.4. The controversy concerning the role of 


pressure fluctuatiors, and momentum transport continues. 


Another conplication woich= 6159 “antroduceds? sf one 
considers momentum transport due to pressure fluctuations is 
that there should be different roughness length scales for 
momentum, heat and water vapour. Beutsaert "19 /S)" nas 
developed equations for the roughness lengths for heat and 
water vapour if the momentum roughness length is known. The 
theory developed is useful for Re w<0. 12 and Rez, where Re 
is the interfacial Reynolds number Re = “AFe Interfacial 


Reynolds numbers found over a water surface include the 


range 0.13<Re <2 where no useful theory exists. 


Since there is no solid theory to use as a guide when 
estimating the roughness lengths for mass, moisture and 
heat, it was assumed that these were equal. Furthermore, 
pressure transport of momentum is assumed to be negligible. 
these Stassumptions “are "consistent?) with” =the ouik-transter 


coefficient formulations of Deardorff. FOr Gonsiscency; 
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Experimental estimates of the bulk-transfer coefficient 
for evaporation have been obtained by various methods. Hage 
(1974) has outlined a number of estimates (Table 3), and a 
wide variation iS apparent among the values for the transfer 
coefficient. Estimates of C Range =fcrom 1 Comm: =. Hage 
argues that the value of 0.0016 for C. , the Lake Hefner 
result, vs vprobably the'tmost  realistac -for Lake ) Wabamun,; 
Since it would best represent the unstable conditions 


present over the lake surface. 


Another route to obtaining a value for the buik- 
transfer coefficient as "to antegrate the profiles of 
temperature, moisture, and momentum based on the Similarity 


theory of Monin and Obukhov, and others. Via this route one 


may ancCcoduce a dependence on Wind, and thermal 
Sula rCation, into the bulk-transfer coefficient 
formulation. By defining faekbulbk-otuanstcr cocificient in 


this manner one should be able to apply the result to any 
wind and stability condition, and obtain meaningful results 


over the heated plumes also. 


The following pages contain the steps by which a 
formulation of the bulk-transfer coefficient based on 


profile theory is obtained (Deardorir W363) 
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Some Estimates for the Bulk-Transfer Coefficient for 
Evaporation (after Hage 1974) 
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Source 


Jacobs (1942,1943) 
Pravet (1960) 
Budyko, et al. 
Wust (Sverdrup, 1951) 
Bunker (1960) 

Deacon and Webb (1962) 


POMC, etn al. Choma) 


Marciano and 
Harbeck (1954) 


Linacre et al. (1970) 


Bveks, 1971 
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CrenG, < rye a result of previous assumptions. 
Accordingly, by integrating (3.3) one obtains: 
3S l 
he Re = | 
u-u a = 3.9 
i lk ae a ae 
5 ° 
Integrating (3.5) one obtains: 
U ee a es 
(0-0, Sere 2 (Sea10} 
(“P5Ee) : 


Again assuming a common roughness for momentun, heat and 


moisture. 
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Similarly, for water vapour in the neutral limit: 

Thus, for this formulation assumed by Deardarft (1968),, 
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also. 


322.3 Extending the Theory to Non-Neutral Stability. 


The stability of the atmosphere is described by the 
Ratsuo ba/Lis LES one Ptconbines = (3276)-Sto 8 (38) "and (3.2) one 
gets: 
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3.2.3.1 Unstable Case 


Assuming that hy OS, in the unstable case, ie., that 
both potential temperature, and specific humidity are 
conservative in the absence of molecular diffusion, and 


using the profile fcxcmulation of Businger (1966): 
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the results: 
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(Ce te ” Res 
n= “VY e 
gale soo otabvemcase 
im the Gase of a stable stratification, Deardorft “has 
used the log-linear profile (Lumley and Panofsky, 1964, 
pi17) to represent the momentum, temperature and humidity 
profiles. 
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Tt. has heen assumed that the neutral transfer coefficients 


aLepegudl-mebguarions (S.72)icand® (3223) yield armuadtacic ain 
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Deardorff has solved this guadratic, and his curves for 
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ts Al and @ =20. Deardorff explains that the 
% 

discontinuities in=the slope ofmthe curves at (Ri), =9 arise 

from the fact that @ = & and 6 = & , and comments 


that more accurate measurements and theory are needed at 


Cause value Of Ris 


Dyer (1975)) suggeste et nar p = 0, = , =145 (z/L) 
and thus the he would all be equal to 5. it is) noted 
that the choice 6 =5 would give a better match for the 
slopes of the Wey), Cubye; under neutrad conditions: 
Purthermore, the results of Businger et ai. (7871) would 
indicate that WA VSP ieardky unity under strong 
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Stapiiicy. thus 8 and oy should be approximately equal. 
The formulas above are cumbersome to use, and because 
the constants in them are not known to any great accuracy, 


an approximate formula given by “Deardorff was usedetar 


computations. In the stable region these formulas are: 
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FIGURE 5. Ratios of the transfer coefficient to its neutral 
value trom diabatic prdotie theory, over) a Wide Tangesof 
Bulk-Richardson numbers (after Deardorff 1968). 
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the bulk-transfer coefficient at neutral stability. 


3.2.4 The Drag-Coefficient at the Limit of Neutral Stability 


In order to apply the previously-developed theory, it 
Wemecleer (that an estimate Pf£or (C, ie (Cu, OL (c,) ie 
required. Near the neutral Ta mact tur bulence-induced 
variations in temperature become small and thus difficult to 
measure accurately. The wind fluctuations, however, remain 
measurable irrespective of the stratification, and thus Ceo). 


Should be easier to determine than(c, ). 


There are -a number of recent estimates of the drag- 
coefficient at the near-neutral limit of stability (Table 3 
aide) Fig: 6). For wind speeds in excess of approximately 4 
Rispethe ftornulas for C, Of (Charnock 9164955, | Smith | and 
Bankes §(1975), Hacks (1972), ‘Shepherd etvad. (A972) ands 
(1969) are in fairly close agreement. The equations are of 
the following forms: 


(Charnock's equation) 
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where the constant C has been estimated as; 
C=11.0 Charnock (1955) 
c=10. 4 Wu (1969) 
C= 1053 Hac Ws ena 72) 


C=10. 6 Smith and Banke (1975) 


) vo . ae 
+ 
‘ 7 - ui 
r . ‘ie te ante 4 
y. | he Bai | 
Deh Ns. ee abs ‘hs hi, 


ak oa 


by Wi} 7 


j jin , 4 sulle 2 oh a st el ; ty dneapinn Perr wad ‘i 
p — a q Ve id ms e 4 na ih 
ba, 4 nighate vet aes isa) \¥ ia ies OF, Sem 


i ut 
¢ wo 45 7 x a i a 
10 ge i Oy a ‘wWeeg“iiwgs we stead 
, p ae : f : 


Pr ) sody Boe bieie per { 1g +6199 0 si ah pe 
vi yh mero ttrorSeudds ba thes ogre \Lovatuone ul 7) 
Ey, . <n 
4 : add 
‘esti iseeaa nid it: Hvedvequsase vide boas 
ai 
' AD : ‘i i ' : | wy i ; 


U H . % si ¢ a 
\ ned ef PTR Fel OF tTesens 1 


4 


eye’ } ‘i! i moteurs: > » TR) ¥ Ay ee a aun fi a Pet -C 


be t ~ = I, _ 
= a ae , a Pe ‘ c's 7 
Oy wen tal bob ved, ‘de: eee ant te Nar 
F Then 


don Bo eas ik ri) jak Qhoayy jaik aot pacha): 


} oe | 
‘ie! : 7 § , 2 ‘Oren ' p 7 v :y or Boy foarat, hs e a 


, bcos? Nae Whe | 
; al”  ,nieeeepe Coote (ole? me 
ih H 3 we P 
at oe ~* = «4 ‘< 4 ie ! i a + oy a : - a) 
| feet? pet, a Le} @ 
vs ; i * 7 
i ' 7 m4 “Sh, ou me 
Pe A 
7 i : é aon 
= AA a pols ae aa ve 
pie po ntas 


a uy rove | | 


PIGURE 6. Drag coefficient and wind speed.. 


(1) Smith and Banke (1975): 
(8.30) mashed curve 
(B23 1) eeline 
(2) Brocks and Krugermeyer (1970) 
(3) Hasse (1968) 
(4) Miller (1964) 
(5) Mayakelet @ilsui1972) 


(6) Shepherd et al (1972) 


(after Smith and Banke 1975) 
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103C) = 0.63 + 0.066a,, (seed) 


Smith and Banke (1975) 


103c, [= 0.36 +40.10u,, (3.30) 


Shepherd et al. (1872) 


Charnock's equation was preferred to the others because it 
was based ona physical argument, and was readily 
transformable from one height to another. This equation has 
only one constant reguiring determination. The Charnock 
equation was chosen to represent the drag-coefficient at 
Wind espeeds over 4 m/s, with a value of 911 “for the 
undetermined constant. This value appears to produce a mean 
functional value for Cy between some of the observed values 


Of iahbtes cand 4 and tFignhe. 


At very low wind speeds, however, one would think that 
the roughness length should approach “that fOr 
aerodynamically-smooth turbulent flow. The measurements of 
Hicks (1972) have indicated that this may be the «case. at. 
remains uncertain as to how the drag-coefficient should 
behaves tingithe: ftransition | reguon from smocth~to-rough 
biurbuLent, ifhow. Indeed, the flow over the waves may be 


transitional at all windspeeds, and the arguments of 
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Charnock may not .be applicable there (Garrat and Hicks, 


1973, Kitaigorodskia, 1968 and others) « 


In this thesis “the  drag-coetficient ts assumed to 
remain continuous in transition from smooth-to-rough 
turbulent flow. At low wind speeds, the flow is represented 


by the profile equation: 
Ug k — (5-33) 
And at higher wind speeds by: 


ROG + Ul Q 
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ix Wee eae 


The drag-coefficient is the greater of that given by (3.33) 
or (3.34) depending upon the wind speed. The variation of 


the drag-coefficient with wind speed appears in Fig. /7. 


ine (82033) and HS. 54) a hot Chew drag=cocrticaent, mad tame s 
usally assumed that the surface drift velocity u, .1s 
negligible, or at least very small when compared to the wind 
Speed Us dicks (1972) measured the surtace dmift, and Lound 
that u, = Use This result was incorporated into the 
calculations to follow. 

Wu (1975) has also examined the surface drift and has 
Obtained values for u equal to 73etom4s or autor etetches 


such as. those of Lake Wabamun. There is little difference 
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O.-O04u, and the error in choosing one over the other is but 


-2% Of u which would result in a 1% error in estimating Cp- 


3.2.6 The Method of Computation 


The equations developed thus far are applicable to 
winds, temperatures and specific humidities measured at one 
level, and water temperature measured at the surface. 
Usually, however, winds and temperatures are not available 
at the same level, as is the case in this experiment. The 
winds were heasured at 4 m above the water surface, and the 
temperature and relative humidity at 1 m to 1.5 m above the 
water. The water temperature was measured at a few 
centimeters below the water surface. This water temperature 
may be somewhat different than the "skin" temperature of the 
water. In the absence of solar radiation, a quantitative 
estimate has been made by Saunders (1967) of the temperature 
difference between "bucket" and actual surface temperatures. 
In this work, no correction was made for the difference 
between the "skin" ce surface temperature and the 
temperatures measured Some cm below the surface. The effect 
of solar radiation on the "skin" temperature is not clearly 
understood, and long-term measurements of radiation are not 
available. Furthermore, the Ryan thermograph was iocated 
near the shore, and, qualitatively, one would expect 
nocturnal water temperatures to be slightly lower than mid- 
lake temperatures, and daytime water temperatures to be 


slightly higher than mid-lake water temperatures. 
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Qualitatively then, the effects of site and depth of water 
temperature measurement would Seem to cancel, and were thus 
ignored for calculation. The magnitude of the "skin-bucket" 
temperature difference is estimated by Saunders (1967) to be 

023 0C in the absence of radiation, and. thus is a 


relatively small quantity. 


In this study, the anemometer and thermohygrograph were 
located at different heights. Provided one uses~ the 
previously-mentioned profiles, one may obtain the following 


system of equations for parameters measured at levels 1 and 
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Ze 
5 ue Claw, ter BLa> te 
(U,-4,) = Ax 
ke Ud - VX vy or Wek ©) 
z, t (2538) 
UU, = 
or awe (2589) 
Thus, assuming that q.°T, «4 es are obtained from 


measurements, one is left with five equations in five 
unknowns u,,U*,U /L and (CS) ae Because an analytic 
solution to these equations did not appear possible, the 
system was solved numerically by successive approximations 
The subroutine for the solution of these equations appears 
in Appendix A. Hourly values of Cor, abet, « (CQ), and the 


evaporation rate were obtained. 
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The method for evaporation calculation described 
previously is rather complicated and laborious. It would be 
preferable to use daily, weekly or monthly averages of 
temperature, wind and humidity to compute evaporation. For 
this reason the evaporation was recomputed using averages. 
Hage (1975) has examined the problem GE computing 
evaporation using monthly mean values of temperature, 
relative humidity and wind in the bulk-transfer equation for 
evaporation. Hage expanded the saturation vapour pressure 
in a Taylor series using the Clausius-Clapeyron equation, 
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relative humidity and wind speed to estimate covariance and 
non-linearity errors. During the summer months it was found 
that the model predicted that these errors would be less 
than about 5%. Hage found that averages could not be used 
during the winter months because wind and temperature were 
correlated because of the passage of synoptic-scale systems. 


An expansion like Hage's, of the stability-dependent 


Ke 


formulation for ~bulk-transter equation was not possible 
because the buik-transfer coefficient is not monotonic with 
respect to wind speed. An estimate was made, however, of 
evaporation using averaged values of wind, temperature and 


humidity and a constant bulk-transfer coefficient. 


The Lake Hefner experiments resulted in a mean valve 
for the bulk-transfer coefficient of 0.0016. It was also 
possible to obtain a constant or weighted-average value for 
the stability-dependent bulk-transfer coefficient over a 


period of time. 


Weighted averages are defined in the following manner: 


Fr 
ta 
f Ce (ul= de 8 Gqts Vaedat - J 
= a ato to eB 
i =u.) C a rAt 
(usu) C= 75) a \w * i cs 
) et (3.40) 
¢, 


Such an average bulk-transfer coefficient may be defined for 


P . = 7 5 ey ee ; n = 
daily, monthly or larger time periods. For this study such 
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a mean was computed over a period of one day, and over all 


the days of measurement. 


3.5 A 


The water-budget method was used to check ° the 
consistency of the results of evaporation calculations. ‘The 


water budget for Lake Wabamun may be written: 


ait) 


Cees sate ett Oe. 0 s (3.41) 


Z e R, G iss 


The runoff was estimated using Thornthwaite procedures 
described by Laycock (1967). The Thornthwaite procedure can 
not be expected to give good results if precipitation is 
distributed very unevenly over the month. Runoff was 
equated to surplus water by the following estimate of 


Storage distribution for the Lake Wabamun basins 


a5 = (eS + Usnloet U4 me Oe eet oS (s540) 
a y, 1 4 © 1o 


Where S, represents surplus water at a storage capacity of 
x inches. The equation above is due to Professor A.H. 
Laycock, Department of Geography, University of Alberta, as 
are the monthly surplus estimates on which this calculation 


was based. 


The net inflow of ground water as well as the net 
evaporation are obtained as residuals in the water-budget 
calculation. The groundwater inflow and outflow are assumed 


to be negligible for the purposes of calculating 
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evaporation. A large groundwater inflow would appear as an 
ever increasing discrepancy between evaporation calculated 
by the water-budget method and evaporation calculated by the 
mMaSS-transfer equations. Laycock (1975) suggested that 
groundwater flow into Lake Wabamun is probably small 


compared with other terms in the water-budget. 


Limits on the amount of enhanced evaporation due to the 
thermal discharge were obtained by considering two extreme 
cases of thermal effiuent behavior. In the first instance 
it was eeeancd that the heated water was instantaneously 
mixed with the waters of the lake, producing a rise in 
overall lake temperature. In the second instance it was 
assumed that no mixing of the plume and lake water occured, 
and that the total excess energy of the plume was dissipated 
by vadiation and turbulent transport of sensible and iatent 
heat. Enhanced evaporation based on variable bulk-transfer 


coefficients and the previous considerations was calculated 


by the procedures outlined below. 
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3.20.2 Two 8 


To model the behaviour of the thermal effluent, the 
lake was assumed initially to be in equilibrium with its 
Surroundings under steady-state conditions. Los vas Jalsa 
assumed that. the heat flux, whether through conduction or 
infiltration of water through the bottom of “the Lake, is 


negligible. The energy budget may then be written as: 


We : Oe ‘ 2h oe (3.43) 

where: 
Bm ey ren Aer AO end (3.44) 
The Of: and Q,, denote solar and terrestrial radiation, 


respectively, and the-radiation is incoming or outgoing as 
indicated» by the jarrowssioThesoutgoing terrestrial, radiation 
QOrseor”, but since € =.97 one approximates spe Tabs 
further calculations. With this approximation {3.44) may be 
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0, = Pa uece OO GRY (pce) 
(3.47) 


ie t 
where CC. > se “tio. (4) (3.48) 


If the rate at which heat is added to the water is Q per 
unit surface area, then the energy budget may be rewritten 


asis 
/ f / 
OP % OP eee (3.49) 


Wirth such an~"added input of heat the water surface 
temperature would rise by an amount 4T, - Assuming foc the 
sake of this approximation that the temperature and specific 
humidity of the air, and the wind speed remain unchanged 
from unperturbed conditions. over this surface, then one may 


Weite: 


t 
OMe=> Ospilsitiy Caphirn Qanekce © cel (3.50) 
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where Or. cape Te AU ce) (3.51) 
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( 
@ = woe Ge C, ¢ sty Wy) WA) 
H (395) 
a ‘ea = 2 3.5! 
an eS OS Cm @) (3.54) 
Combining the equations above one obtains: 
DO = SIG Bese a GA a 
( ‘ 
+ fo Sy OOD |< are Set 
( 
nerd Cp rtp Cate rh ae citady seats “ee 
(SR SS) 


For the case of the thermal effluent diffusing throughout 
the whole lake, a further simplification is possible. Since 
4 ft is very small compared to Toe and Cp On one can expand 
the Spectiiceghinadity sandr eradialion “0Sihnge fir st=order 
approximations with respect to temperature. By doing so one 


obtains: 
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{4 6 fs {au (3.56) 


ve es 100 


4 
ie C |<. + Ve 4, . 5B9716*S4%87.6 Jonn seen 


(3-127) 
where E,= input of thermal energy and A = lake area. 


The equation above may be used to obtain an estimate of 
the rise in lake temperature for a wuniformly-dispersed 
thermal effluent. The relative magnitudes of the terms give 
the partitioning of energy between O, + Get and 0 Thus, 
for a given undisturbed lake temperature, it is possible to 
determine the enhanced evaporation knowing the energy input 


tobthe lakevand ithecnean wind. 


When one considers the case of a plume of uniforn 
surface temperature Lab T, it, eGuti pram | with wits 
Surroundings then, neglecting heat-_cansfer Copthe west of 
the lake, the plume size becomes a function of the wind 
speed, provided T., Te R,, OT, and E are known. Once again 
Hees possible to obtain e partitioning “OL Vener gy between 
the latent, sensible, and tLadiative components of energy 


transfer from the plume surface. 
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4.1 Introduction 


In this chapter the eresults of the different methods of 
calculating evaporation from Lake Wabamun are presented and 
discussed. Before such results can he presented, however, 
ib is ‘mecessary to consider the |/relzability of the 
measurements. Table 4 gives the estimated errors of 


measurement of the field instrumentation. 


There are major sources of errors other than the 
intruments. Marciano and Harbeck (1954) found that mid-lake 
winds were 10% higher than winds measured at the shoreline. 
By mounting the wind-measuring instruments on bermS or 
protuberances extending into the lake, and then using 
records from stations for which the wind was off the lake, 
it was hoped that measurements so obtained would be more 
representative of mid-lake conditions. Because Lake Wabamun 
Vemeenetaliy: warner ‘(Wan the Soverlying sat, thers e215 
probably a thermally-driven convergent component in the wind 
field around the lake. Thus, even by taking wind 
measurements downwind from the lake, there exists the 
possibility that the winds will be underestimates of the 
true mid-lake wind. Practical considerations dictated 


mounting wind-measuring instruments on top of berms or piers 
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extending into the lake. Unfortunately, the physical 
Presence Or such @asnound) on theitlat surface’ of the water 
would produce an increase in wind speed over the mound, and 
temperatures and specific humidities measured there would be 
representatives Ore Melqntse, Slightly (lower: @thanmsactual iy 
measured (Fig 8). An estimate of the greatest error which a 
ridge the size of the Sundance berm could produce is 
obtainable from the inviscid-flow approximation for flow 
perpendicular to a cylinder. From these considerations it 
aS-estimated that the winds could be as much as 2% higher 
over the Sundance berm than those measured at anemometer 
level away from the berm. The anemometers were calibrated 
in awind tunnel, but over the lake the flow was turbulent, 
and one would expect a degree of overspeeding of the cup 
anemometers (MacCready 1966). No compensation was made in 
Thee Calculations: Or caja bration Oe ycorrec. nace & such 
overspeeding. The magnitude of such overspeeding would 
depend on the amount of turbulence present, and hence would 
be affected to some extent by the stability of the flow. 


The magnitude of Such an overestimate is expected to he less 


than 10%, however. 


The air temperature and humidity aS measured by the 
thermohygrograph mounted in a Stevenson Screen may also be 
in error due to screen-lag in adjusting to the environment. 
Bryant (1968) found that the lag time decreased from 30 
minutes under calm conditions to 6.5 minutes with a 7 m/s 


wind. In this study no compensation was made for screen 
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tag= Errors due to the neglect ~of ‘screen lag may have 
Ltesulted in systematic differences in the calibration of the 
Sundance and Wabamun thermohygrographs. The weekly 
calibration measurements were taken shortly before noon at 
Sundance, and in the mid or late afternoon at Wabamun for 
most of the checks made of these stations during the summer 


NOntchSs in) 19745 


Probably the greatest source of error in the water- 
budget calculations is the question of whether or not the 
rain-gauge measurements are representative of the lake as a 
whole. During the summer of 1974 radar observations at the 
Edmonton International Airport confirmed that sShover or 
thundershower activity occurred on all days for which there 
were precipitation and radar observations available. It is 
noted that radar pictures were available only when the radar 
was working, and when the meteorological technician was not 
EOOSE DUS eWEth @OThCre Gdulves se forethese ) Gascon (ne badd t 
data were rather incomplete. For water-budget calculations 
an average of the precipitation measurements from the two 
stations is used to estimate average lake precipitation. A 
Dazred-t© "test applied! to "the® #8 ™day “total precipitation 
produced confidence limits of 16.4 43.1 cm at the 90% 
confidence level (16.4 46.3 cm at 95% level), assuming that 
the errors in measurement were randomly distributed. The 
distribution of rainfall may not be even over the lake due 
to the interaction of the wind and terrain. There may also 


be systematic errors in the catch of the tipping-bucket 
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rain-gauges but the magnitudes of these remain undetermined. 


The modified Thornthwaithe procedure (Laycock 1967) 
based on monthly PACA MARAIRIOR and temperature data failed 
to indicate any runoff from the Lake Wabamun drainage basin. 
The validity of this procedure when applied to monthly 
amounts is guestionable during two periods of heavy rainfall 
{cay 10) Vandy day 930). When resuits from water-budget 
estimates of evaporation are examined there appears to he 
ENNOLt BLOU up to, 10 days atiter the heaviest raantatl.. 9 8 
period of runoff was suggested by a period of 10 days during 
which the totaled evaporation decreased. Another possible 
Source of error in water-budget calculations was set-up of 
themiake suctace. | A COrLection was) included in lake -level 


Calculations to compensate for set-up of the water surfac 


ra) 


(Tppen 1966). It was found that the set-up was rather smaitli 
most of the time. The effect of waves and Seiches was not 
determined, although it is thought that seiches would induce 
reading errors in the estimation of the mean water level 
from the chart recording. Waves could induce systematic 
errors in water level if the Stevens water-level recorders 


did not respond linearly to variations in water depth. 


4.2.1 The Bulk-Transfer Coefficients 


The bulk-transfer coefficients and evaporation were 


ie 1 
A ‘ uM Ap 2 le ie i) 
ase ame Teale, meat g tittaae ag 
i S : a 
’ \W a “ re i iy 7 iy 44 nat aia ray, sf 
if pond oe wnat Henley yee 
sb i. ve if x é | 


> ay Tae Hen “aaptie’ ory “4 etds | Lom ¥: 


i i | chia s it ry? F ae iy + a 
Ana ’ or wnt bah * fdaoe Pai laa’ pe Z ad vce, 
7 , i 2, ve 7) nes, 1 _ q 
“ a f ; le nl 4 r - 
ison ie! ine. sh Ae 
: Wot aw ay .' ‘= 


ao ei ; a 
7 Dok > 
we ie FS, a i wn! A a 
ae. reheat eolas: Feehan: ite ry ot i Leneah mip: 


f ‘ : ae y 26 a io 
Lye “a Tten t+ DPPRLANSS. Was pip IM Ave, m4 wo aes 
i op ? 7 - 


mr ae pete 
nei ned ¢ yoo Ta ryeh Of es ci aot: ae 


; ; rawr ae Cen met eL they vos ay 
he Solo Bb id Setesogua 2391 ou nate ne 2m 
wm 7 L : td 7 ‘ it} : 
‘is i ' we Te : = | a 1s aan 28 > 
f ee cal t of! b es ‘5 ht 1% . 4 tv u bso AJ 3 , » iy ds | sls P if | [ 


ha 
i 


‘, o7 4 Ae a _ 
. ws p. f : f a b ie a s 1 7 re it 8 7 » A 
{ j 4 ; Lae ; are 2 Fe By . rn 1 / £ i‘ aD ana tet did 
: ; re | 
, : LG) 
‘ 7. ~poe tal seemiffieta, ad. 
a. Toe F ’ ” ~ ae x 
i ‘a * ) , ; .: ast 
“ge ea Paow .e 7 rT [ iy Phi } ‘a i ‘7 ed 


Ten! ry 
4 i a sls > 

i ne pe) te i + ge 
bat Bove Bae Shweta Het eta oad 
mr n> 


4 “sss Ps ae ay. oo Ay 
Loe, et? pi fz épvode ta 
‘ ra u a ; robe 
ee eel ree Be ete ee 
(coe Ott) Pee alone: oi). “a RAC 
7 : i 1) ie hats - 


«DRS BS RST 2 
S| nk 


55 


Calculated from data reduced to hourly values using the 
theory and procedures outlined in Chapter ERED The 
coefficients, and evaporation were then chosen from either 
or both stations depending upon whether or not the wind was 
TLORe = the=-watcreatesthesstation.-+ Theshouphyevalueseuot the 
bulk-transfer coefficient at thermograph height are plotted 
in Fig. 9. A considerable diurnal variation was found for 
the value of the bulk-transfer coefficient. Weighted daily 
averages of the bulk-transfer coefficient for evaporation 
appear in «Fig. 10% A comparison of these average 
coefficients with the daily mean air-water temperature 
difference (Fig. 11) shows a rough correlation as expected. 
A mass-transfer weighted average value of the bulk-transfer 
coefficient for evaporation of 0.00162 was computed for the 
78 days of the experiment. This result appears to be in 
close agreement with the results of Marciano and Harbeck 
‘Gis is It must be remembered, however, that the average 
bulk-transfer coefficient which was calculated is valid ata 
height of approximately 1.2 m, and that, strictly speaking, 
the @Hetner formula piss valid for winds and Specific 
humidities measured at 10 m. An attempt was made to compute 
an average 10m bulk-transfer coefficient using the weight 
function outlined in Chapter IIf. A weight funetaon  ~ for 
See was obtained by taking the ratio of evaporation to the 
computed 10 m bulk-transfer coefficient. Unfortunately the 
approximate formula used to calculate the buik-transfer 


coefficient vroduces an unrealistic weight factor Ons 
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Figure 9. Hourly selected C, at 1,2 m. 
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peELodst yofAestablies Pstratittd cations This was because for 
Ri>1, both evaporation and the bulk-transfer coefficient for 
evaporation were small. However, the ratio of these two 
near-zero quantities produces a weight function which is 
many orders of magnitude larger than weight factors during 
lapse conditions. It suffices to say that one should not 
try to divide zero by zero. Weight functions at thermograph 
height were not plagued by this problem because the wind was 
interpolated to this height, and the weight function there 
was obtained by the product of wind and air—-water specific 
humidity difference. The approximation was better at this 
height because (Rif was usually (less than) 1. The 
evaporation formula that results in a bulk-transfter 


CocumrTctent 070200162 45: 


(= = Cc. pP Cu-ad lq ah (4.1) 


wheras the Hefner formula is: 


ee ca Ga Oy 43) (4.2) 


Because u, = 0.035 u the average bulk-transfer coefficient 
derived from Lake Wabamun data becomes 0.00156 in terms of 
the Lake Hefner equation. Most drag-coefficient 
determinations, including the ones from which the neutral- 
stability drag coefficient was obtained for this thesis, 
e.g-, Smith and Banke (1975), have obtained the drag 
coefficient over water without compensating for surface 


drift. By compensating for surtacah dcitt-, “whe terditral— 
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Stability drag coefficients so determined would have been 
higher by approximately 8%. thiis= sand the Lack for the 
Surface drift in the Lake Hefner formulation account for 
about 12% of the observed difference in the evaporation 


estimated by the two formulas. 


Hage (1974)) has asedy 0.0016 as the value for the bul k— 
transfer coefficient For evaporation to calculate 
evaporation throughout the open-water Season. Such an 
appLoach snould bee Gustrived@ as long ~as\ thew sans liey 
remains close to that observed in the Lake Hefner experiment 
at other times in the year. The temperature record Shows 
unstable conditions to exist during the ice-free period and 
air-water temperature differences like those observed during 


the period of this study. 


A cunning total of the evaporation was computed as a 
Change in the Jake level using the stebility—-and-wind— 


dependent formulation of the bulk-transfer coefficients. 


Mromacesults, are plotted Sin) fag. aes Because water 
temperature data were unavailable because of recorder 
malfunctions, there are two breaks in this record. The 


evaporation and bulk-transfer coefficients were chosen from 
either observing station or an average of both, depending 
upon whether or not the wind was blowing off the lake at the 


observing site. This was done in an effort to ensure that 
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the instruments and measurements were consistently well 
Within the boundary “layer of the lake, ‘and thus® the 


conditions for a constant-flux layer would be fulfilled. 


Evaporation was also calculated using the same data and 
the Lake-Hefner equation with a constant bulk-transfer 
coOstefrcrent of, 020016, (Figs) 13,714 -and /15))-) Because vind 
speed was not measured at 10% m aS is) required sforueth ies 
fomulation, the wind at 10 m was extrapolated using the wind 
profile formulation of the previous chapter. For unstable 
cOnpdrTevons the ~Ccorrection » from 4s 9m) to 10m ease asmaly 
(2 BOE u,)- However, LOL Stable conditions the 
correction can be large. Under stable conditions the wind- 
profile formulation is not particularly good, because the 
assumptrton that Z/L be siall as not sfuliilied. | The specitic 
humidity was not corrected to 10 m in this caiculation. 
This would iead to a Slight underestimate of the air-water 
specific humidity difference and the evaporation. 
Evaporation calculated using the constant | bulk-transfer 
‘coefficient is plotted as total evaporation ten dS- ss oeeetO 
15. Again, gaps in evaporation calculations are due to 
missing water temperature data. From these graphs of total 
evaporation it is evident that there is little difference in 
using either Wabamun or Sundance observing station, or a 


choice from those stations depending upon whether or not the 


wind is from the lake. 


Some bias is possible in all of the mass-transfer 
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estimates used because of the shoreline location of the 
thermohygrographs. One would expect that warmer water would 
tend to remain near the surface, and as a result the warmest 
water would tend to be dragged by the wind and advected to 
the lee side of the lake. There was some evidence for the 
occurrence of this from temperature profile measurements 


taken oy Nuttall (1974) This effect, “and the pessipigity os 


calibration errors, would tend £0 Paccoun wm sore etne 
differences in evaporation calculated at Sundance and 
Wabamun. If the wind tends to be onsShore at Sundance a 


cneaes portion of the time, this station would presumably 
be within the water-vapour blanket of the lake for a larger 
portion of time. The specific humidity at Sundance would be 
higher and the surface specific humidity would be Jowered at 
Wabamun due to upwelling of cooler water. fis WOuld 


partially explain any differences between the two stations. 


Rather than uSe hourly values, it would be more 
convenient to use daily, weekly or monthly mean winds, 
humidities and temperatures to compute lake evaporation. 
The wind-and-stability-dependent formulation is clearly 
unsuitable for estimating evaporation uSing mean winds. 
This formulation, or at least the wind-dependent part of the 


formulation, is not monotonic with wind speed. Daily mean 


wind speeds generally would fall on the part of the drag- 
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coetticient curve  — (Fig. 7) which is representative of 
smooth turbulent flow. The wind-and-stability-dependent 
formulation was not used for estimating evaporation using 


means of atmospheric variables. 


Because ce obtained from the stability-and-wind- 
dependent formulation approximated the 0.0016 value given by 
Marciano and Harbeck so closely, .-this constant value was 
used to estimate evaporation from Lake Wabamun uSing mean 
values of air and water temperatures (Fig 16). It was noted 
that the use of mean data produces a 10% underestimate of 
evaporation compared to the use of hourly values. 
Evaporation was also calculated using the "min-max" method 
Of sSpring . and Schaefer -(1973),% butwthis method Mroduced 
rather low values (Fig. 17). This method used the Hefner 
equation, and the daily minimmum RH and maximum temperature 


to compute the daily air-water specific humidity difference 


and the evaporation. 


Estimates Based on the Water Budget Method 


——— a 
——- oe ae ee ee 


omparison with Other Estimates 


—-— ——— == 


The water budget was used to check the accuracy of the 
preceding methods of calculating evaporation, and to infer 
whether or not there may have been a net infiow or outflow 
of groundwaters to or from Lake Wabamun. To begin with, it 
was assumed that seepage and runoff were negligible. 


According to the modified Thornthwaithe procedure of Chapter 
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itt no runoff should have occurred during the course of the 
experiment. The remaining quantities in the water-budget 
equation are changes in the water level, the precipitation 
and the discharge through Wabamun Creek. These quantities 
are plotted in Figs. 18,19 and 20, respectively. The 
evaporation is plotted as totalized evaporation because the 
reading C€rrors of the Chart recorder tend to canveel) cach 
Otheryon Such a plot.” it appears that after the major vain 
Onpeday 30) there “was a period of runoff of about) 10 days 
duration. At other times the assumption that runoff was 


negligible appears reasonable. 


The spikes on the plot of total evaporation as a 
funcewoue Of etl mea(Fig.s §24) “aremditficult stom explain. BY 
Gomparing the total evaporation (Fig. 21) with the total 
Ratmidtte (hig. 19) ait appears that these spikes are 
Correlated with rainfall events. A LTikely explanations 
that strong winds associated with these showers produced 
seiches which made reading of lake level difficult, and 
introduced noise into the graph of totalized evaporation. 
This is corroborated by the wind speed observed either prior 
to, or during the time when Such spikes appeared on the 


water-budget estimate of evaporation. 


Using procedures outlined in Chapter III, an estimate 
for a mean bulk-transfer coefficient was calculated from 
water-budget evaporation, the air-water specific humidity 


difference, and the wind speed. Only the periods from day 1 


a) " a q oo i” 
Ho bomaeede sen he 


Lie 


| age Plans 
yan ok ww, Mt as J 
gig vs = " ae are. = 7 
i Smee ( 


| 4 r pr eat Lov sa fa eA pind it 
; 7 aie a 
: ¥ - eine * i goaab 
- Aoum) tm ut, 4 Ware aia oan 
: , i‘. <7 i) ry ch _) lL id 
| mw Gh’ ibihas’ . ‘a ae faa as “ft na? 
Rio oe ee it 
Ur (Area me ee aha Oe | 
noven Seuttecot Bs Bal t0ly ef aot teow 
yon wkene ety ho. anna oe ‘eashaye 
f : ry 
49% +o 19 ae they oy dh bs «2924 e dows “NO as ny 
af | WER 
t | ¢ ehyeq 2 rémv wrod? oe cnh 1 
: i” . ; ; ch 1” a W AJ 
+ ug mis sek | An’ ore ae i o dA ef 4028 = 7 b 
i ee , \ < ‘ 
Aula 
ehidenesge% ateceue ext Rett pea 
= ; ea Cae 
é ; 7 wy ! : 
oe min fa nee) A *) Aq va? he rer iz i ait 
"4 r 7 *- id vey, t 
- TA ty 
fe US POR anit Yo yok spa 
mats? 
ae 05 i%) Ke Loane ge? satot wiht _ ph (36 
7 


AAT BR TA ij F rs (er 
yionas seteeve LLatateg as ty 
2 = tivo 4 aa 25 7 pene - ot ky i 
i oy fy i. ”, ® a 2 ae, oe ‘ ap * 0 


! : oe) eu a i 
townl $464): gh Loa e223 oben . 
ny 
‘ a : 


nee ye. ; 
3 Pach 


% Wgyie: nS 
. ys ) 


ee an ead 
a. i 


~) i 


: ae 7 ate A 
190A Pa, qe dupe ne 
wr ea Mee 


nu See 
eth 


iit 


*TOAVT 29eeM UNWeGeM ONeyT “ST ein3TtyY 
shep Ut uty 


08 02 O09 OS OF DE 02g 


| RR hc 


Ot 


COSSoL 


STS SE h 


(CISY W) [ese 1907eM 


Vato 


SE ae GL 


i zn 


os 


\ aan en aly a | 


As 
F ‘ : ; e. F 
“ : 5 " - 1 
* 1 7 ™“‘. e ¢ 
} : a - q 


72 


‘uoTyejatdtoaerd j[e}O}R unuleqemM sHeT “ST ean3ty 
sep ul duty 


08 OL Og OS OF OF 06 OF C 
SS SS ee eee 


00"0 


¥3°0 


26° 0 


+... 


Ov°O 


(uw) uotjzeitdtaerg 


lie) 


*yderig UnuleqemM YsnNOAUI esBzeYOSTpP VATIeTNUND *OZ7 V<AN3BTY 
sAep UT WT] 


08 OL 09 OS OF O€ 0d OT 0 
ee ee ee 


§0°0 70° 0 


Past 
(W) ad1eYISTp dATJeTNUIND 


St au 


} 
1 
ws 8 RA ie RSS 


02°0 


Oh ae ie 
if eo a : i i Qn 


' } 


74 


*poyzeul Zespng-197eM JU WOAZ peje[NoTeo UoTJeATOdeAD sATReTNUND °TzZ eANn8Ty 
: shep UT owty 


08 OL 09 OS OY GE 02 OT 
a aes tlle ate a Nn 


g0°G 


SH 


vo" 


Crs 


(uw) uotzet0odeaq 


» 
ba Ss 
ag aft 
4 


ee 


a 
* 
» 


ny 


pe £7CEL<fF 
2 
. 


, Pe ; 
bes . , a : 
= Woe 
rf / ~ S - ata ata 4 f 3 » “we io! pana wet 
4 ft ial . : 
ws! » Te 1 Fi a * ; rm 
a . 4 ¢ @2ea ls Ue 
” | eer 
b é | ie 
ha > ea) 
Oy & *' = 7 - =) 4 ‘ee 
id 4 i i 7 \ % , i 
; J oe am 4 
a ( a ; a Po = ; A Bis 
; a 
oe : ma, cae 
‘ : i 
- 


¥ 


¥ 
~ 


| anne 
TIS 40 
2 
ie ae 
4 
¥s 


; 
Py 7 : : a ay Maa 
: ve 7 = 7 ri 
= : a : ae 
f iy 
b : Whe! 
> 
ia - rs \ ‘" 4 os ‘ ~ ; ah 
\ & Y a : 
ae 


a2 


dha 


to 17 and from day 40 to 78 were used in the calculations 
because there was a possibility of runoff ro davyeovtosday 
40 and there was a period of missing water temperature ‘data 
arters? dayen 7. Not much confidence is felt for the early 
part of the water-budget evaporation estimate, since there 
waS another period where a runoff could possibly have 
occured aroundidaysl0c%1etThe bulk-transter’ *c0effieient ’ § for 
evaporation was estimated at 0.0016+0.0006 for days 1 to 17, 
and OC 016f0. 0004S for Fdays GO%EOr 738.9 Ther error linits on 
the former figures refer to 90% confidence limits on 
rainfall obtained from a paired-t test statistic on rainfall 
at Sundance and Wabamun. In addition to this, the water- 
budget estimate should have an additional uncertainty of at 
least 20% since the bulk-transfer coefficient is calculated 
uSing the integrated product of wind and air-water specific 
humidity difference. The evaporation and mean bulk-transfer 
coefficients computed from the water-budget are about the 
same as the evaporation and the mean bulk-transfer 
coefficients of the mass-transfer methods used. Runoff from 
themesnallvacreeks®Sarctnd’ the® lake 7oeret “Arneluded) “in tie 
eaicalataons,* “vould “serve to make the water-budget 


evaporation and bulk-transfer coefficients slightly greater. 


Considering the confidence limits which the 
precipitation places on the evaporation, the differences 
between mass-transfer estimates of evaporation, and water- 
budget estimates of evaporation, At is not ‘possible * to’ say 


confidently whether there LS at weet anftiow -orvoutriow or 
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Water fromethe Lake sbottom.. <Previous studies «(Nursall, ..et 
ali, 19 725 gander Est zeand akrouse, 6197/3), have also estimated 
evaporation from Lake Wabamun. WNursall et al. (1972) used 
monthly means to calculate evaporation from Lake Wabamun and 
obtained runoff using the method of Bruce and Clark (1966). 
BaucesandsClark estimate, runoff, as. .f%. 08 precipitation 
falling on the catchment area of the lake. Hage (1974) 
found this formula to underestimate snowmelt runoff and to 
overestimate runoff from showery summer rainfall. A 
modified Thornthwaite procedure (Laycock, 1967) was used at 
estimate runoff from Lake Wabamun. The method predicted no 
runoff; but, as»was suggested by Laycock. (1975), there could 
still be runoff associated with a Single large rainfall. 
Laycock also noted that the modified Thornthwaithe procedure 
could be computed for shorter time periods, and if this were 
done  wehdate it, could ppredict.srunofttiseassoccUrired aiter the 
Storm, rainfall on day 30., However, because there ,»appeared 
tom Demanpossibdiity-ofsrunofft only after, thassonc event, sot 
was decided to neglect runoff completely in the water-budget 
calculations, and to concentrate on those times when it was 


reasonably certain that there was no runoff. 


Previous estimates of evaporation from Lake Wabamun 
vary as to estimates of yearly evaporation amounts, e.g, 
Wursali, et al. (1972) 51.96 cm, Hage (1974) 54.6 cm, and 
Ferguson (1973) 69 cm. The first two represent evaporation 
estimates based on monthly averages of atmospheric variables 


and the bulk-transfer method, while the latter uses the 
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Meyer formula, and an unknown data base. The preceding 
results of this chapter indicate that the use of mean values 
for computing evaporation produces an underestimate of 
evaporation, when compared to the evaporation computed using 
hourly values (10% underestimate). This underestimate could 
be due to the different covariances between winds and 
temperatures at land and water sites. The evaporation 
figures of Ferguson appear to be rather high in comparison 
to the other mass-transfer estimates of evaporation. Hage 
(1974) used monthly mean winds from the Stony Plain 
observing site, about 32 km east of the lake. Hage thought 
this station to be more closely representative of monthly 
wind speeds at Wabamun. During the course of this study, 
however, it was found that, over the period of 78 days, the 
unadjusted mean wind speedS at Wabamun and Sundance were 
very close to the mean wind speed. at the Edmonton 
miteinataonal “Airport, observing site, (25C80 in/s, a2. co oem s 
and 2.869 m/s, respectively). Hage (1974) used the monthly 
mean wind speed ratios Edmonton International: Edmonton Stony 
Plain, and Edmonton Industrial:Edmonton Stony Plain, to 
adjust Edmonton International and Edmonton Industrial winds 
to those at Edmonton Stony Plain when the latter were 
unavailable (prior to 1966). Hage (1974) found that wind 
speeds at Stony Plain were lower than those at Edmonton, and 
assumed these would be more representative of the Wabamun 
area. The results quoted here would tend to question this 


assumption. However, it must be remembered that this study 
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coveredsa perwod fofonly pasie days, Sand that, -thas .sumper 


season may not be representative of the climatic normal. 


The possibility of a net outflow of groundwater from 
the lake remains. This possibility was to some extent 
Bein LOrced Dye thow .1ndings Of NUursallveet eal md 1972 \erand 
REeatzeand) —KbLouse, (197305 Because of the isotepic 
composition of the lake, these authors conclude that the 
lake is under a strong evaporitic regime. However, 
estimates retenie evaporation, precipitation, sand .nuneLe 
computed by these authors require a net intlow of 
groundwater into the lake. In order to explain the low 
salinity of Lake Wabamun, these authors postulate that there 
must also be a discharge of groundwater from the lake. The 
Study “of Fritz and’ “‘Krouse (1973)) indicates’ “that wells 
drilled to the south-west of the lake contain water for 
which the isotopic composition 1S similar to the lake water. 
Based on geological evidence (Carlson, 1970) they conclude 
that a buried Pleistocene riverbed acts aS a recharge and 
discharge aquifer for Lake Wabamun. They postulate that the 
aquifer Links Lake Wabamun with Lake Isle to the north-west, 
is supplying Lake Wabamun with water of the right isotopic 
composition and lower salinity, and is flushing the lake 
through groundwater to keep the lake at a low level of 
salinity. The mass-transfer evaporation estimates made in 
this thesis, when compared to water-hudget evaporation 
estimates, indicate that there is little if any net flow of 


. is he eee 
groundwater either to or from the lake. Furthermore, it is 
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possible that Lake Wabamun is flushed sufficiently by 
periodic outflows of water from Wabamun Creek , for which 
discharge data were incorporated into the water-budget 
Calculations. The creek iS intermittent in nature, and 
discharges water to the North Saskatchewan River whenever 
the water level on Lake Wabamun is sufficiently high. Hage 
(1974) found that discharge of lake waters ocurred for lake 
devetse urgher= than, 722-/ mo. The findings wouldpanutearce 
that seepage, if it exists, is taking water out of Lake 
Wabamun and that this seepage, along with periodic 
discharges of lake water from Wabamun creek, account for the 
PSOLOPIC *cOmpositton and Salinity of the lake; The fact 
that there iS an intermittent discharge from Lake Wabamun 
indicates that, over longer time periods, the lake treceives 
more precipitation and runoff than is evaporated or seeps 
Phrough the Lake bottom. Assuming a “normal” climate, -cne 
could equate this discharge to the allowable evaporation 
enhancement due to power development. 


4.5 Evaporation Enhancement due to Thermal Effluent 


=e —— _ — 


Lake evaporation can be thought of as a mechanism by 
which the lake gives up energy absorbed as solar radiation. 
Neglecting the exceptional case where energy is advected to 
the lake, the incident absorbed solar energy is in 
equilibrium or near equilibrium with the long-wave 
terrestrial radiation emission from the lake, and the latent 


and sensible heat given up to the atmosphere over the long 
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term. As was formulated in Chapter III any additional 
energy input to the lake is also partitioned among these 
components of the energy budget. The power plants on Lake 
Wabamun put a thermal load of 270 MW into the lake when 
Gpereting vatetheit gmaximug capacity of 900 MW.) | ifeall this 
power, were to be, used exclusively for evaporating water 
there would be an increase of lake evaporation of 4.1 cm per 
year. This figure is considerably lower than the 5.8 to 8.1 
cm maximum possible yearly evaporation calculated | by 
Ferguson (1973) using the Meyer formula. The high figure 
may be due in part to an overestimate of wintertime 
evaporation kased on the assumption of too large an open 
water area by the power plants. This figure may aiso be 
affected by Ferguson's interpolation of evaporation data 
from evaporating-pan stations not representative of the Lake 


Wabamun area 


A comment is made at this point that it is conceivable 
that the open-water area may not be due wholly to the 
increased heat input of the power plants, but may be a 
result of the mechanical mixing of the warmer water below 
the 1ce by the kinetic energy of (he water—jev “from the 
plants to produce an open-water area. This effect would 
take out heat during the winter months at the expense of 


energy stored during the open-water Season. 


Hage (1974) has estimated Winter and summer ilake 


evaporation on a monthly basis, and has obtained a yearly 
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evaporation enhancement Of gee Lae ls Hage used crude 
estimates of open water and plume areas from the work of 
Hogg (1973), Nursall and Gallup (1971), and Nuttall (1974). 
The LANDSAT satellite pictures (Figs. .22 and 23) Show the 
wintertime open-water alLedas “at “different. times) in the 
Winter. Nuttal (1974) has. on a few occasions measured water 
temperatures within the open-water areas. He found that the 
plume occupied only a small portion of the open-water area. 
However, these results may not be representative since there 
was considerable wind on the few occasions that these 
temperature measurements were made. One would imagine that 
the stress exerted by the wind on the water would serve to 
mix the plume and lake water, resulting in a more compact 
plume. Lal and Rajaratnan (1975b) have determined 
experimentally the size of a thermal plume in a_ tank as a 
function of a plume Richardson Number. The results were of 
Interest pe eDUtIy.Ol dd.) sNOte sagLeecraVver yews lent Che tne few 
measurements of the size of the plume taken by Hogg (1973). 
Not knowing what the appropriate water temperature 
distribution should be, it was decided to consider two 
extreme cases of thermal effluent behaviour as outiined in 
Chapter IIi, and to compute the relative magnitudes of 
enhanced energy transfer by radiation and latent and 
sensible heat transfers. Because the energy input is known, 
this analysis should set upper and lower limits on 
evaporation enhancement. Figs. 24 ta-26 show the felative 


percentages of the energy transfer components, and the rise 
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FIGURE 22. LANDSAT satellite pictures of Lake Wabamun. 
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satellite pictures of Lake Wabamun. 
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in temperature of Lake Wabamun as a function of wad speed. 
Pie we) look= ct a typical @wind §speed of & w/s, theslake 
teomperaturesis Seen fo rise, Thy 0218 CG) for an “unpolluted 
surface temperature of 0 C, and .12 C in summer for a water- 
surface temperature of 20 C. The percentage of extra heat 
input dissipated by latent heat rises from 30% in winter to 
58% in summer. A crude, all-season average gives a figure 
of 44% for the portion of total energy input by power plants 


that is dissipated as latent heat. 


Figs. 27 to 29 show the relative percentages of the 
energy transfers, and the size of the thermal plume required 
to dissipate the excess energy of the power plants, assuming 
a plume temperature 10 C higher than the water temperature. 
In this case, latent heat accounts for 35% of energy loss in 
Winter and 60% in summer. A crude, all~Season average gives 
48% as the percentage of energy input that is dissipated as 


latent heat. 


From the results above one would expect that latent 
Heat accounts for slightly lessee chan, S040 scl US thes theanat 
pollution from the power plants, and that enhanced 


evaporation shouid be no more than about 2 cm/year. 


Because Lake Wabamun occasionally discharges water from 
Wabamun Creek during years with high water levels, it should 
be possible to estimate an average discharge of water from 
the lake. If the normal discharge is greater than that 


required to alter the lake level by 2 cm/year, then there 
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should be little or no problem of enhanced evaporation 
associated with the operationwot the power plants at.) thesr 


present capacity. 
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CHAPTERRY 


CONCLUSIONS 


The period of time covered by this study was short, and 
restricted to the summer months. However, it is doubtful 
that a longer period of study would yield a efinement of 
the conclusion that the wind-and-stability-dependent version 
Oteethegebulk—-transtergecoefEicient edoessnottyield eapbebter 
evaporation estimate than a constant bulk-transfer 
coefficient. Within the uncertainties imposed by measured 
components of the water budget, either estimate produces 


equally valid results. 


It was encouraging to find that the time-averaged value 
for the bulk-transfer coefficient determined from profile 
theory was close to the 0.0016 value from the Lake Hefner 
studies which was used by Hage (1974) in his estimation of 
Lake Wabamun evaporation. it is also concluded that the use 
of daily averages of temperature and humidity data from 
Fdmonton International Airport was justified to within the 
gecuracy of thismrstudy. While these results were 
encouraging, it is necessary to point out that much of the 
theory used to construct the wind-and-stability-dependent 
bulk-transfer coefficient requires some validation and 
refinement. In particular, the scaling lengths z, and az, , 


y 


and their interdependence, need to be known explicitly over 
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a water surface. The stability dependence of R75 needseto 


be determined in a more accurate and unambiguous manner. 


A limit for enhanced evaporation due to the injection 
of heated waters was*cdlculated. &t was found that a little 
less than 50% of the thermal load on the lake was used to 
evaporate water, The remaining 50% is dissipated as sensible 
neat and) Terrestrials  fadiation. With a power-generating 
Capacity of 900 MW, lake evaporation is increased by at most 
2 cm yearly. A mean discharge of 6 cm annually was computed 
for Wabamun Creek. Taking this discharge to be surplus to 
that reguired to maintain the present lake level, it should 
be possible to double or triple the present power generating 
Capacity without having to import water to maintain lake 
levels. Lake~level decreases would be greater during dry 
years. However, the decreases could be at most double the 
present lake-level decreases. Perhaps, for this reason, and 
because of the impacts of the hot plume on the local plant 
and animal life, the added expenses incurred by Calgary 
power in diverting North Saskatchewan River water to cooling 


ponds is reascnable and justified. 
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LIst OF SYMBOLS 


DESCRIPTION 


area 

bulk-transfer coefficient, drag 
coefficient, and bulk-transfer 
coefficients for heat and water vapour, 
respectively 


transfer coefficient under neutral 
conditions 


Specific heat of ale tac OnsS can 
pressure 


evaporation 

energy input to thermal plume 
acceleration due to gravity 
precipitation, runoffiand ground. water 
inflow expressed in terms of lake-level 


change, respectively 


eddy diffusivity for momentum, heat and 
water vapour, respectively 


von Karman's constant 

Monin-Obukhov Length 

evaporation from the lake surface, 
subsurface outflow and streamflow fron 
the lake, respectively 

atmospheric pressure 

total heats flowgapheat. Eicwe duei, to 
cadiative, latent and sensible heat 


transfers, respectively 


Sola Gc and terrestrial radiative 


‘transfers 


specific humidity, saturation specific 
humidity of a water surface, fluctuating 
and at a level z, respecively 
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/ 7 / 
U,Vy,W,U VV gw 
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DESCRIPTION 


Reynolds and interfacial 
numbers, respectively 


Reynolds 


Bulk Richardson number 

relative humidity 

time 

temperature, surface, fluctuating, 
Virtual, fluctuating virtual and 


Cemperatune pat dslevel zienesvecruively 


wind at the surface andeat usa meleve lume 


respectively 


horizontal and vertical wind components 
and fluctuations of these, respectively 


friction velocity 


height above ground (m) 


roughness, or scaling lengths LOG 
momentum, moisture and heat, 
respectively 

emisivity 


change in surface elevation of the lake 
Kinematic viscosity 
air density 


Stefan-Boltzman constant 


shear stress and surface Shem AS 
respectively 

dimensionless gradients _ of wind, 
temperature and water vapour, 


respectively 
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List of FORTRAN Variables 


Variable Description 

TTH the height above water of the 
themohygrograrh 

OMQS air-water specific humidity difference 

cv virtual temperature 

UAN wind speed at anemometer height 

ZTH height of the thermohygrograph above 
water 

ANHYDF difference in heights between the 
anemometer and thermohygrogtcaph 

CEOVCH the ratio Ce/C,, (=1) 

CDN the neutral-stability drag coefficient 

RI the bulk-Richardson number 

CD the drag coefficient 

CE the buik-transfer coefficient for 
evapocation 

UTH wind speed at thermohygrograph height 
relative to the water surface 

TH the water surface temperature 

EVAP the evaporation rate expressed as a lake 


level change in m/s 
u10 the wind speed at 10 a 


CE10 the bulk-transfer coefficient for 
evaporation at 10 an 


OBML the Monin~-Obukhov length 
CDN10 the neutral-stability 10 m drag 
coefficient 


CD10 ™ the 10° “drag coefficient 
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